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Net-baryon number fluctuations in (2+l)-flavor QCD 
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We present a lattice study of net-baryon number fluctuations in (2+l)-flavor QCD. The 
results are based on a Taylor expansion of the pressure with respect to the baryon chemical 
potential. We calculate higher moments of the net-baryon number fluctuations and compare 
with the corresponding resonance gas results. We find that for temperature below 0.9T C 
the fluctuations seem to agree with the hadron resonance gas predictions. Close to T c , 
higher moments are increasingly more sensitive to the critical behavior of the QCD phase 
transition. Furthermore, we estimate the radius of convergence of the Taylor series as well 
as the curvature of the transition line in the temperature chemical potential plane. 

§1. Introduction 



Fluctuations of conserved charges, as baryon number and strangeness are gen- 
erally considered to be sensitive indicators for the structure of the thermal medium 
that is produced in heavy ion collisions!^ 1 In fact, if at non-vanishing baryon number 
a critical point exists in the QCD phase diagram, this will be signaled by divergent 
fluctuations of all quantities that can be connected to the fluctuations of the chiral 
■ | order parameter as, e.g., the baryon number density!^ ^ 

We present here results from lattice calculations of net-baryon number fluctua- 

-C tions in QCD with dynamical light and strange quark degrees of freedom. The results 

are based on calculations with an improved staggered fermion action (p4-action) that 
strongly reduces lattice cutoff effects at high temperature and are performed as joined 
work of the RBC-Bielefeld Collaboration. 4 ' The values of the quark masses used in 

\^q | this calculation are almost physical; the strange quark mass, m s , is fixed to its phys- 

t— I ■ ical value while the light up and down quark masses are taken to be degenerate and 

equal to m s /W. This corresponds to a pion mass of m^ ~ 220 MeV. The lattice 
spacing of the lattices we analyze here is with a ~ 0.25 fm (N T = 4) rather large and 
we expect a cutoff dependence due to a deformed hadronic spectrum. 5 ' However, for 
the quantities we discuss, which are ratios of moments of net-baryon number fluctua- 
tions, the information on details of the spectrum cancels in HRG model calculations. 

"^ ■ It thus may be expected that lattice cutoff effects arising from a distortion of the 

spectrum are less severe in such observables. They may be directly comparable to 
the hadron resonance gas results and to heavy ion experiments. 

§2. The Taylor expansion method 

Direct lattice calculations at nonzero baryon density are impossible by means of 
standard Monte Carlo methods. We follow here the Taylor expansion approach as 
described in detail in Rev. [6]) ,[7]). Starting from an expansion of the logarithm of the 

typeset using VTpJ^.ds (Ver.0.9) 
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QCD partition function, i.e. the pressure, one obtains 



where 



rpA 


~ VT3 


n z. 


\V,±,HB) 


i 


.-¥) 


cf(T) = 


-x B = 

i\ Xl 


1 

i! 


Qi 


In Z(V,T, 


Vb) 


d(j*B/T)* 


VT 3 





(2-1) 



(2-2) 

Mb=0 

Here [Ib is the baryon chemical potential, which can be obtained as appropriate 
linear combination of the quark chemical potentials. Note that we treat the up 
and down quarks as degenerated, both in mass as well as net-quark number den- 
sity. Accordingly, isospin and electric charge chemical potential vanish. Due to 
charge conjugation symmetry only even coefficients are nonzero. The Taylor expan- 
sion coefficients cf are calculated on the lattice as expectation values at [ib = 0. 
Derivatives of the partition function with respect to \xb are also know as moments 
of net-baryon number fluctuations. The second derivative defines the quadratic fluc- 
tuation X2 — 2!c^, while higher derivatives give higher moments (xf, Xe 1 • • • )• 

§3. The hadron resonance gas 

Within the HRG model interactions are encoded in the thermal creation of 
hadronic resonances, in accordance with their Boltzmann factor. 8 ' The HRG is 
known to describe the observed particle abundances in heavy ion collisions very 
successfully® and has recently also shown to describe particle fluctuations observed 
at RHICP® For the heavy baryons we can utilize Boltzmann approximation and 
obtain for their contribution to pressure 

% = % (y) 2 K2{mi/T)coMB ^ B,T) ■ (3 ' 1} 

Here di, mi are isospin degeneracy factor and mass of baryon species i, respectively. 
From Eq. (j31j) we can easily see that the /^-dependence factorizes and that ratios 
of moments of baryon number fluctuations are either unity {xf /X2 = X% /X2 = 1) 
or are simply given by their ^-dependence, e.g. we have X2 /X\ = coth(/is/T), 
X f/x 2 B = tanh(^/T).m) 

§4. Results on net-baryon number fluctuations 

In Fig. mieft) we show the 4 and 6 order net-baryon number fluctuations, xf 
and Xbi normalized by the quadratic fluctuations \ B - F° r temperatures T<0.9T C 
the fluctuations seem to agree with the predicted HRG result which is unity for both 
quantities. For temperatures T>0.9T C the fluctuations start to deviate from the 
HRG model. As expected x B increase more rapidly close to T c , since this quantity 
diverges in the chiral limit, where as Xa wm develop a kinkA^ In Fig. [T^right) we 
show the quadratic net-baryon number fluctuations at [ib > 0. We compare the 4 th 
order (data points) and 6 th order (dashed lines) results. We find, that errors from 
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Fig. 1. Shown are 4 th and 6 th order net-baryon number fluctuations, normalized by the quadratic 
fluctuations (left) and quadratic fluctuations at [jlb/T > (right) as a function of temperature. 



the truncation of the series (|2-ip become sizable for pb/T ~ 1.5. Finally, in Fig. [2] 
we plot our results as function of the center of mass energy y/s along the chemical 
freeze-out curve by using a suitable parameterization!^ These results can be di- 




Fig. 2. Shown are ratios of moments of net-baryon number fluctuations along the freeze-out curve, 
as a function of the center of mass energy ^fs. 

rectly compared to the net-proton number fluctuation that are measured in heavy 
ion experiments. Note that the ratios have been expanded to obtain their leading 
order ^-dependence. We expect that for center of mass energies of ■ v /i<20 GeV, 
i.e. /Ub/T> 1.5, higher orders become significant. It is thus important for the under- 
standing of future heavy ion experiments to obtain Xs to good accuracy from first 
principal lattice calculations. 

§5. The radius of convergence 



It has been argued that the radius of convergence can be used to determine 
the location of the QCD critical point Jffl-LLSJ'LIJJ 1 One way to define the radius of 
convergence of the pressure series (p[p/T 4 ]) (Eq. ()2-ip ) is by 



lim p n with p n = pp '/T = Jc^/c^ +2 . 

n-^-oo V 



(5-1) 



The radius of convergence can be estimated in a similar manner by the coefficients 
in the series of the quadratic net-baryon number fluctuations (/o[x^D- Each order 
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Fig. 3. Estimates of the radius of convergence as explained in the text (left), same data plotted in 
the (T, /i)-diagram, together with the freeze-out curve and a band for the transition line (right). 

p n will than differ by a constant factor which approaches unity in the limit n — > oo. 
In Fig. |3](left) we plot the first two orders of both estimators. We find, that at 
T = 0.95T C the estimated radius of convergence is consistent with what we find from 
analyzing the truncation errors, i.e. we find p$ ~ (1.5— 2.5)p,g/T. In Fig. fright) we 
plot the same data in the (T, ^s)-plane and estimate the curvature of the critical line 
from p^, assuming that the radius of convergence is limited by the phase transition. 
Interestingly the lower limit of the curvature which we estimate from p^ is roughly 
in agreement with what was found from an analysis of the O(N) universal scaling!^* 
Note that all the results that have been presented here need to be confirmed on 
finer lattices. First iV T = 6 results on \2 an d xf have already been published. 7 ^ 
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